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Summary 

This report describes the development of an interaction data base and a numerical solution 
to the transport of baryons through an arbitrary shield material based on a straight ahead 
approximation of the Boltzmann equation. The code is most accurate for continuous-energy 
boundary values but gives reasonable results for discrete spectra at the boundary using even 
a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H 2 O). The 
resulting computer code is self-contained, efficient, and easy to use. The code requires only a 
very small fraction of the computer resources required for Monte Carlo codes. 

1. Introduction 

As NASA continues to develop a vigorous space program, tools for the analysis of optimum 
shielding against space radiation are a continuing requirement. The tools must ultimately 
account for the very complex mixture of radiations in the space environment and be complete 
in the physical description of the processes involved. An incomplete model must, by necessity, 
have restricted capabilities that are not always appreciated at the engineering level and may 
cause errors in vehicle design. Still, a complete model must be computationally efficient in order 
to provide a useful tool for design work. In order for the model to be used with confidence, 
some effort toward model validation must be made, _ 

Space contains the most complicated mixture of diverse radiation components known. When 
these components interact with materials, many new and varied radiations are produced. This 
places enormous demands on tools for model development analysis. Furthermore, a common 
basis for assessing risk from such an environment is in itself a challenge to model development. 
The model must likewise allow for all potential elemental materials and allow inhomogeneous 
configurations. The present work is the beginning of this task. 

Monte Carlo computer codes have been written that meet many of the above requirements 
(ref. 1). However, the enormous computational requirements have caused their use to be avoided 
in the space program. In an earlier report (ref. 2), we presented a relatively complete transport 
code for high-energy nucleons. The data base for that code was complete but somewhat 
inaccurate. The purpose of the present report is to describe both the improved computer 
programs developed for the calculation of the transport and the interaction of high-energy 
nucleons (baryons) with materials. The methods, based on the direct solution of the Boltzmann 
equation, have been developed over the last several years (refs. 2 to 5). The present goal is 
to document a relatively complete description of the basic physical processes and an improved 
input data base. Future work will concentrate on improving the data base and validating the 
computational procedures. 

2. Theoretical Considerations — The Boltzmann Equation 

In moving through bulk material, particles give up energy to the medium through 
atomic/molecular and nuclear interaction. These processes are described by a Boltzmann- 
like equation that we use in a time-independent form. The equations in the straight ahead 
approximation to be solved (ref. 1) are 
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where <f>j(x,E) is the differential flux density of type j particle at x with energy E; S(E) is 
the proton stopping power; ut is the ion-range scaling parameter; cr p (E) and cr n (E) are proton 



and neutron total cross sections, respectively; and E f ) represents the differential cross 

sections for elastic and nonelastic processes. It is convenient to define new field quantities as 
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which can be rewritten as integral equations with the boundary at x = 0. The following results 
are given: 
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The functions and coefficients of equations (2.1) to (2.12) are presented in the next section. 

3. Transport Coefficients 


3.1. Stopping Power 

In passing through a material, an ion loses the larger fraction of its energy to electronic 
excitation of the material. Although a satisfactory theory of high-energy, ion-electron inter- 
action is available in the form of Bethe’s theory utilizing the Born approximation, an equally 
satisfactory theory for low energies is not available. Bethe’s high-energy approximation to the 
energy loss per unit path (that is, stopping power) is given as 
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where Z p is the projectile charge, N is the number of target molecules per unit volume, Zt is 
the number of electrons per target molecule, m is the electron mass, v is the projectile velocity, 
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/3 = v/c, c is the velocity of light, C is the velocity-dependent shell correction term (ref. 6), 
and It is the mean excitation energy given by 

Z t ln(J*) = £/„ln(£„) (3.1.2) 

n 

where f n represents the electric dipole oscillator strengths of the target and E n represents the 
corresponding excitation energies. Note that the sum in equation (3.1.2) includes discrete and 
continuum levels. Empirically, it has been observed that molecular stopping power is reasonably 
approximated by the sum of the corresponding empirically derived “atomic” stopping powers 
for which equations (3.1.1) and (3.1.2) imply that 


Z ln(7) = 

3 


(3.1.3) 


where Z and / pertain to the molecule, Zj and Ij are the corresponding atomic values, and 
rij represents the stoichiometric coefficients. This additive rule (eq. (3.1.3)) is usually called 
Bragg’s rule (ref. 7). 

Sources of deviations from Bragg’s additive rule for molecules and for the condensed phase 
are discussed by Platzman (ref. 8). Aside from shifts in excitation energies and adjustments 
in line strengths as a result of molecular bonding, new terms in the stopping power appear to 
be due to coupling between vibrational and rotational modes. Additionally, in the condensed 
phase, some discrete transitions are moved into the continuum, and collective modes among 
valence electrons in adjacent atoms produce new terms in the absorption spectrum that must 
be considered. It is usually assumed that the experimentally observed additive rule shows that 
molecular stopping power is the sum of atomic processes. In contrast, Platzman proposes that 
molecular bond shifts for covalent bonded molecules are relatively independent of the molecular 
combination, thus resulting in an additive rule. On the basis of such arguments, Platzman 
suggested that ionic bonded substances should be studied as a rigid test of the additive rule 
because of the radical difference in bonding type. He further estimated that ionic bond shifts 
could change the stopping power by as much as 50 percent. Recent results of molecular bond 
shifts on mean excitation energies are discussed in references 9 to 11. Effects of the physical 
state have likewise been studied (ref. 12). 

The electronic stopping power for protons is adequately described by equation (3.1.1) 
for energies above 500 keV for which the shell or “tight binding” correction C makes an 
important contribution below 10 MeV (ref. 13). For proton energies below 500 keV, charge 
exchange (electron transfer) reactions alter the proton charge over much of its path, so that 
equation (3.1.1) is to be understood in terms of an average over the proton charge states. 
Normally, an average over the charge states is introduced into equation (3.1.1) so that the 
effective charge is the root-mean-square ion charge and not the average ion charge. At any ion 
energy, charge equilibrium is established very quickly in all materials. Utilizing the effective 
charge in equation (3.1.1) appears to make only a modest improvement below 500 keV, an 
indication presumably of the failure of this theory based on an empirical basis (refs. 13 and 
14). We have utilized empirical fits to the proton data; the resultant stopping power for protons 
in water is shown in comparison with the evaluated data to Bichsel (ref. 15) in figure 1. 

The electronic stopping power for alpha particles requires terms in equation (3.1.1) of higher 
order in the projectile charge Z p resulting from corrections to the Born approximation. The 
alpha stopping power cannot be related to the proton stopping power through the effective 
charges. Parametric fits to experimental data are given by Ziegler in reference 16 for all elements 
in both the gaseous and condensed phases. 

The electronic stopping power for heavier ions is related to the alpha stopping power through 
the corresponding effective charges. The effective charge suggested by Barkas (ref. 17) is used: 
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The total stopping power Sj is obtained by summing the electronic and nuclear contribu- 
tions. Other processes of energy transfer such as Bremsstrahlung and pair production are 
unimportant. 

For energies above a few MeV/nucleon, Bethe’s equation is adequate provided that appro- 
priate corrections to Bragg’s rule (refs. 9 to 11), shell corrections (refs. 6, 13, and 14), and 
an effective charge are included. Electronic stopping power for protons is calculated from the 
parametric formulas of Andersen and Ziegler (ref. 13). The calculated stopping power for pro- 
tons above a few MeV in water is shown in figure 1 in comparison with data given by Bichsel 
(ref. 15). 

Because alpha stopping power is not derivable from the proton stopping power formula using 
the effective charge at low energy, the parametric fits to empirical alpha stopping powers given 
by Ziegler (ref. 16) are used. Applying his results for condensed-phase water poorly represented 
the data of references 19 and 20. Considering that physical state and molecular binding effects 
are most important for hydrogen (ref. 9), the water stopping power was approximated by 
using the condensed-phase parameters for hydrogen and the gas-phase parameters for oxygen 
(which are known experimentally). These results are compared with experimental data for 
condensed-phase water (refs. 19 and 20) in figure 1. It appears that Ziegler overestimated the 
condensed- phase effects for oxygen, since the gas-phase oxygen data gave satisfactory results 
as seen in figure 1. 

Electronic stopping powers for ions with a charge greater than 2 are related to the alpha 
stopping power through the effective charge given by equation (3.1.4). For water, the condensed- 
phase formula of Ziegler for alpha particles gives probably the best stopping powers for heavier 
ions. Calculated results for 1 -’0 and 56 Fe ions in water are shown in figure 1 in comparison 
with the Northcliffe and Schilling data (ref. 21). Good agreement with Northeliffe and Schilling 
for 56 Fe ions is especially important since their data seem to agree with the range experiments 
of J. H. Chan in Lexan material, a polycarbonate resin developed by The General Electric 
Company (ref. 22). The stopping powers in Lexan resin and tissue-equivalent material can be 
calculated in a way similar to the procedure given above for calculating the stopping powers in 
water. 
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Although the above methods are the best yet available, they do not adequately represent 
the data at all energies and elements. The energy loss per unit path length is given by 

= 4 g ( 3 . 1 . 8 ) 

mr 

where an approximate representation of true stopping number B is given in the braces of 
equation (3.1.1). The determination of S e then rests upon the accurate knowledge of the mean 
excitation energy I t and the shell corrections C. In practice, one invokes some sort of parameter 
fitting involving the experimental data on stopping power and the quantities It and C. 

We have recently initiated a different approach of calculating the stopping power of 
atoms. The main thrust of this approach is to calculate exactly, in the context of the Born 
approximation and nonrelativistic wave functions, the stopping number B with no assumptions 
such as those underlying the Bethe theory leading to equation (3.1.1). Thus, the knowledge of 
B would rest on knowing the radial integrals for the process of excitation as well as ionization 
when a projectile passes through matter. To the best of our knowledge, this approach seems not 
to have been tried even for the helium atom. We have recently calculated the radial integral 
for the optically allowed transitions in the He atom and helium-like ions using the screened 
hydrogenic model. The model describes the atom by single-particle hydrogenic wave functions 
and treats the initial state and the final state by two different, effective charge parameters Z, 
and Zt, respectively. The model is able to reasonably reproduce the existing dipole oscillator 
strength values with little effort, and nonrelativistic numerical values for bound-bound and for 
bound-continuum transitions are available for many target helium-like ions. The model has 
also been successful in reproducing the known dipole polarizability values and in predicting 
the other unknown values. Once the radial integral is evaluated for all momentum transfers, 
it is an easy matter to obtain the stopping power of the helium atom for a projectile such as a 
proton or a heavy ion. The same approach then could be extended to include other atoms. 

This approach is an ambitious undertaking but is more satisfying in that the calculations 
are made directly for each atom, thereby avoiding the inherent approximations underlying the 
Bethe equation (3.1.1). Thus, the calculations do not involve the apparent fittings involving 
parameters such as It and C. Furthermore, the ejected electron distribution in energy and 
angle as well as the atom excitation spectrum should also be available through this approach. 

3.2. Total Nuclear Cross Sections 

After many decades of experimental activity at various accelerators with ever-increasing 
energies, the cross sections for two-nucleon interactions are reasonably well-defined. Although 
recent advances in the theory of the two-nucleon interaction in terms of phenomenological 
meson exchange models (ref. 23) show considerable success, a simple parameterization of the 
experimental data is sufficient for our purposes. For E > 25 MeV , the proton-proton (pp) total 
cross section ( mb ) is found to be reasonably approximated by 

<jpp(E) = ^1 + {40 + 109cos(0.199v^)exp [— 0.451( J E7 - 25)°' 258 ] j (3.2.1) 

and for lower energies, by 

<Tpp(E) = e 6 - 51e_<E/134) ° 7 (3.2.2) 

These forms are compared with experiments (ref. 24) in figure 2. For E > 0.1 MeV, the 
neutron-proton ( np ) cross section is taken as 

a np {E) = 38 + 12 500 exp [-1.187(E - O.l) 035 ] (3.2.3) 

and at lower energies, by 

a np {E) = 26000e~ (£ /°- 282) ° 3 (3.2.4) 
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These forms are compared with experiments (ref. 24) in figure 3. 

The low-energy, neutron-nucleus total cross sections exhibit a complicated fine-resonance 
structure over a broad, slowly varying background. This background is marked by very broad 
Ramsauer resonances that persist even to neutron energies of 100 MeV. Although a simple 
fundamental theory for the Ramsauer resonances is not available, a semiempirical formalism is 
given by Angeli and Csikai (refs. 25 and 26). Their formalism starts with the usual partial- wave 
expansion as 


a t = 2?rA 2 ^(2£ + 1)[1 - Re(%)] 

l 

(3.2.5) 

with 

m = e i6 ‘ 

(3.2.6) 

where Sg is the complex phase shift for the ith partial wave and R e(Z) denotes the real part of 
Z. In the opaque nucleus model, the fact that n^wl for all values of t > Rj A, where R is the 
nuclear radius, leads Angeli and Csikai to assume that 

crt fts 2n(R + A) 2 [l - Re(r/)] 

(3.2.7) 

where rj = 0 gives the usual, opaque nucleus result such that 


R e(r?) = e _ hn(<$) cos[Re(<5)] 


= p cos 

(3.2.8) 

is a reasonable starting point to parameterize the total cross sections, where Im(<5) denotes the 
imaginary part of 5. Their complete parameterization is 

at = 27r(ro^y 3 + A) 2 [a - pcos^yiy 3 — r)] 

(3.2.9) 

where tq = 1.4 fm, and the neutron wavelength is 


4.55 + 1 

/E A t 

(3.2.10) 

The parameters of Angeli and Csikai are adequately approximated by 


1 

" ~ 1 + [2/(3. 8E + 0.1 E^/E + 0.1 E 3 \fE)} 

(3.2.11) 

p = 0.15 - 0.0066 \[E 

(3.2.12) 

q = 2.72 - 0.203 \[E 

(3.2.13) 

r = min{— 5.3 + 1 mVE] 1.3} 

(3.2.14) 


Strictly speaking, equations (3.2.9) to (3.2.14) apply only to At > 40 and 0.5 < E < 40 MeV. A 
simple extension to all values of A% and 0.1 < E < 100 MeV gives qualitatively similar results 
to the experimental data and provides a starting point to representing the total cross section. 
The cross sections given by equations (3.2.9) to (3.2.14) are shown in figure 4. This should be 
compared with the experimental data (ref. 27) shown in figure 5. Note that the data in figure 5 
have only the broad resonances shown. The very narrow resonances have been averaged. We 
now seek some modification to the Angeli- Csikai cross sections to better approximate the total 
cross sections. 
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Our modifications to the Angeli-Csikai formalism are as follows: 

1. If At > 75, then a is taken as 0.18 for values of equation (3.2.11) less than 0.18. 

2. The value of p is taken to be greater than 0.4a unless At > 76 for which p can be as 
small as 0.3a. 

3. A modifying factor of 1 + De~ aE is used with 


f 0.5 (145 < A t < 235) 
\ 1.0 (Otherwise) 


f 1.0 (205 < A t < 235) 
\ 2.0 (Otherwise) 


4. An additional modifying factor is applied as 


Fiji - 0.5 exp [-( A t - 63.54) 2 /2o] 

- 0.45 exp [-( A t - 58.71) 2 /4] exp(-2F) + F 2 } 

where 

F _ f 0.7 {A t < 63; E < 0.8) 

1 \ 1.0 (Otherwise) 

f 0 (E > 0.5) 

F 2 = < -4.95e" 18£ (40 < A t < 42) 
l-1.79e“ 15E (32 < At < 34) 

5. If At < 30, then numerical interpolation between experimental values is used. 

The final cross sections as modified above are shown in figure 6 and should be compared with 
figure 5. 

The total cross sections above 100 MeV have been taken from reference 28. The high-energy 
cross sections of reference 28 have been approximated by 


c r tot {At,E ) = 52.5^ 758 [l + (0.8 + 2.4e -yl ‘/ 30 )e- £ / 135 sin 0 £ ] (3.2.15) 

where the phase angle is given by 

14.41 (E < 40 MeV) ) 

1.29 ln 2 (F) - 7T (E > 40 MeV) / 1 ' 

The expressions (3.2.15) and (3.2.16) are shown in comparison with the theory of Townsend, 
Wilson, et al. (ref. 28) and a compilation of experiments in figures 7 to 10. Equations (3.2.9) to 
(3.2.14) are connected smoothly at 70 MeV to the results of equations (3.2.15) and (3.2.16) at 
130 MeV with an assumed exponential dependence on energy. The total cross section is used 
to calculate the scattering cross section as 



<Js{E) = CTtot(F) - <T abs (E) (3.2.17) 

The total (tot) neutron-nucleus cross section is shown in comparison with experimental data 
(ref. 27) in figures 11 to 14. (Experimental data are shown as the dashed curve.) Also shown 
are the cross sections (listed as “prior”) used in reference 2. Clearly, the present result is a 
great improvement. 

3.3. Nuclear- Absorption Cross Sections 

Qualitatively, the nuclear-absorption cross sections show an energy dependence similar to 
that observed for the total nuclear cross sections. An analytic formula for protons was derived 
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by Letaw et al. (ref. 29) by first fitting the cross sections of Bobchenko et al. (ref. 30) with the 
formula 

a A = 45A? 7 {1 + 0.016sin[5.3 - 2.631n(A t )]} (3.3.1) 

where At is the mass number of the target nucleus. Equation (3.3.1) reproduces the Bobchenko 
data to within ±2 percent (ref. 29). A somewhat better fit to the Bobchenko data is given by 

a A = 45A?- 7 (1 - 0.018 sin 0^) (3.3.2) 


where the angle is 


& A = 2.94 In {A t ) + 0.63sin[3.921n(A*) - 2.329] - 0.176 (3.3.3) 


Equation (3.3.2) fits the Bobchenko data to within the 1.2-percent difference, which is on the 
order of the quoted experimental uncertainty (ref. 30). Although the Bobchenko data represent 
a consistent set of measurements for many different targets and probably define well the 
A-dependence of the high-energy cross sections, they may nonetheless be in error in absolute 
value as suggested by many other independent experiments (refs. 31 and 32). 

Letaw et al (ref. 29) assume the energy dependence for all nuclei to be the same and to be 
approximated by 

f(E) = 1 - 0.62e~ E/20 ° sin(10.9£T°- 28 ) (3.3.4) 

where the nucleon kinetic energy is in units of MeV. We observe oscillations according to the 
quantum-mechanical calculations of Townsend and Wilson (ref. 31) with phase angle 


0 E = 


1.44 (E < 25 MeV) | 

1.331 n(E) - 2.84 (Otherwise) J 


(3.3.5) 


but with an A-dependent amplitude given by 


f{E) = 1(0.3£T°- 22 + O.76e _ ' E/135 )(O.4 + O.9e _y4 ‘/ 3 °)sin0 £ ; (3.3.6) 


The absorption cross section as given by equations (3.3.3), (3.3.5), and (3.3.6) is shown in 
comparison with the fit of Letaw et al. and experimental results in figures 15 to 19. As one can 
see from the figures, it is difficult to assign a figure of merit to the fit, since great scatter in the 
data obscures the result. Generally, above 20 MeV the results are on the order of ±10 percent 
accurate as estimated from the scatter in the experiments. 

Below 20 MeV, the neutron cross sections are represented by numerical data sets at discrete 
energies of 1, 3, 5, 10, 14, and 20 MeV as taken from references 27, 32, and 33. Interpolated 
values between data points at the available target masses are shown in figures 20 to 25. 
Intermediate energy values are found according to 

<7 (At, E) = a(A t , Ei)e-< E ~ E ^ (3.3.7) 


where E, and a are taken according to the appropriate subinterval. The cross sections are 
assumed to be zero at energies below 0.5 MeV. The absorption cross sections for elements from 
lithium to plutonium for energies between 1 and 100 MeV are displayed in figure 26. 


3.4. Fragmentation Cross Sections 

The local distribution of ions and radicals produced in ionizing radiation events is known to 
be an indicator of biological response. The fact that such distributions for high-energy nuclear 
radiation are vastly altered by local nuclear- reaction events has been studied in nuclear emulsion 
(refs. 34 and 35) and is a regular component in risk assessments in high-energy neutron and 
proton radiation fields (refs. 36 and 37). Risk assessments have generally depended on the 
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results of calculational models of these reactions, since the detailed study of such reactions was 
largely inaccessible to experimental study until the advent of high-energy, heavy ion beams. 

The first detailed, relativistic, heavy ion beam experiments were performed by the Heckman 
group at the Lawrence Berkeley Laboratory (LBL) (refs. 38, 39, and 40) in which beams of 
carbon and oxygen were fragmented on a series of targets ranging from hydrogen to lead. The 
momentum distribution of the projectile fragments relative to the projectile rest frame was 
measured for all the isotopes produced. These results will be analyzed to ascertain relevant 
biological factors with their corresponding implications on radiation-risk assessment in high- 
energy, nucleonic radiation fields. An ion fragmentation model will be recommended for use in 
radiological protection and studies. 

Individual nuclear constituents are ejected in the collision of high-energy neutrons and 
protons by direct collision (ref. 41). The remaining nuclear structure is left in an excited state 
that seeks an equilibrium minimum-energy configuration through particle emission (ref. 42). 
This is the basis of Rudstam’s study of the systematics of spallation products produced in 
such collisions in which he assumes that the resultant isotopes are distributed in a bell-shaped 
distribution near the nuclear stability line. The total change in nuclear mass and the dependence 
on the incident projectile energy are treated empirically in Rudstam’s formalism. 

The fragment charge distribution for a given fragment mass Aj is given as 


}{Zf) = exp(pAf —r\Zj - sA f +vA 2 f \) 
where the coefficients show a slight energy and fragment-mass dependence as 

T = 11.8^J 0 - 45 
5 = 0.486 
v = 3.8 x 1(T 4 


f 20 E~ 0 - 77 (E < 2100 Me V) ) 

P ~ \ 0.056 (E > 2100 MeV) J 

where E is the nucleon energy. The complete Rudstam cross section is given by 

<r(A f ,Z f ) = [F^pA^fiZf)] /D 

where 


D = 1.79 


e pA ‘ ( 

03\ 

i V 

pAt) 


O 03' 
M + pAt. 


F\ = 5.18 exp(— 0.25 + 0.0074.4*) 


(3.4.1) 

(3.4.2) 

(3.4.3) 

(3.4.4) 

(3.4.5) 

(3.4.6) 

(3.4.7) 

(3.4.8) 


F 2 


exp(1.73 — 0.0071E) (E < 240 MeV) ) 
1 ( E > 240 MeV) J 


(3.4.9) 


We have applied a simple mass-dependent correction factor to Rudstam’s formula as shown 
in table 3.4.1 and renormalized his cross sections to the total absorption cross section. Many 
corrective factors have been added to Rudstam’s formalism by Silberberg, Tsao, et al. (ref. 43). 
Estimates have also been made by Guzik (ref. 44) for some of the isotopes produced in 
connection with cosmic-ray propagation studies with some attempts at experimental verification 
(ref. 45). 

From a nuclear-model point of view, isotope production at low energy results from the 
formation of a compound nuclear state that decays through particle emission. At higher 
energies, the direct ejection of particles from the nucleus becomes important, and intranuclear 
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cascades represented as sequences of two body scatterings within the nucleus with Pauli blocking 
are the usual means of evaluation (refs. 46 and 47). Subsequent to the cascade, the residual 
nucleus is assumed to be in thermal equilibrium and seeks to minimize its internal energy 
through particle emission (ref. 46). 

The measurement of isotope-production cross sections at proton accelerators does not allow 
the direct observat ion of the fragment products. Customary measurements used 7 or j3 counting 
techniques to identify the isotopes produced. Stable and short-lived isotopes produced in the 
reactions were either not observed or their number was greatly distorted by loss through decay. 
This is particularly true for light-mass targets such as those that are important to biological 
health considerations. Consequently, the fragmentation of carbon and oxygen nuclei by protons 
remained shrouded in experimental obscurity until the advent of heavy ion accelerators. 

One of the earliest experiments performed at the LBL Bevatron, when the ions of carbon and 
oxygen could be accelerated to relativistic energies, used detectors able to measure the energy 
and charge of an ion beam in conjunction with a bending magnet for momentum analysis 
(ref. 38). In this way, the density in phase space was measured for each isotope produced in 
collision with a fixed target. 

The isobar cross sections (olbl) measured by Lindstrom et al. (ref. 40) for 2.1 A GeV oxygen 
fragmentation on hydrogen targets are given in table 3.4.2 in comparison with the results of 
Bertini (ref, 47), Rudstam (ref. 42), and Silberberg, Tsao, et al. (ref. 43). Note that the 
Rudstam results contain the correction factors in table 3.4.1 and are renormalized as described 
above. 

The oxygen-fragmentation cross sections as represented by three parametric forms are 
shown in figures 27 to 31 in comparison with the Bertini results and various experiments. 
The baryon-15 isobaric cross sections in figure 27 show that experiments favor the curve 
of Silberberg, Tsao, et al. Although the Bertini model provides an overestimate, the other 
parametric curves provide improved estimates compared with the Bertini code. The baryon-14 
isobaric experimental cross sections are in reasonable agreement with the three parametric 
curves as well as with the Bertini model as seen in figure 28. Again, the experiments show no 
clear advantage of one parametric curve over another for the baryon-13 cross section as seen in 
figure 29, although the Bertini results appear somewhat low. We show experimental results for 
baryon numbers between 9 and 13 of Lindstrom et al. in table 3.4.2. Clearly, the equally good 
agreement for the Rudstam parameterization and the Silberberg, Tsao, et al. parameterization 
is obtained by baryon numbers 12, 11, and 10. The Bertini cross section is far too low to 
represent the cross section for baryon-11. The baryon-9 cross sections are shown in figure 30. 
(The results of Yiou are reported in ref. 45.) The Silberberg, Tsao, et al. parameterization 
is too high by a factor of 2 or more. The haryon-7 cross sections are shown in figure 31. At 
energies below 300 MeV, the baryon-7 results of Silberberg, Tsao, et al. are favored. 

The measurements of Lindstrom et al. (ref. 40) for relativistic carbon beams are shown 
in comparison with the results from Rudstam (ref. 42) and Silberberg, Tsao, et al. (ref. 43) 
in table 3.4.3 for two beam energies. The good agreement with the results of Silberberg, 
Tsao, et al. is no surprise, since their parameterization was fit to these experimental data sets. 
Note, however, that the Silberberg, Tsao, et al. cross section for mass 8 fragments needs to be 
suppressed. 

3.5. Differential Nuclear Cross Sections 

3,5-1. Nucleon- nucleon spectrum - The nucleon- nucleon differential 
represented (ref. 48) by 

f{E,E') = B ^e-B(E'-E) + e -B£j j ^ 1 _ e-BE'') 

where 

B = 2toc 2 6/10 6 


cross sections are 

(3.5.1) 

(3.5.2) 


I II! 
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In the above equation, me 2 is the nucleon rest energy (938 MeV) and b is the usual slope 
parameter, given by (in units of GeV~ 2 ) 


f 3 + 14e E Z 200 (For pp) 1 
\ 3.5 + 30e~ fi '/200 ( For np j j 


(3.5.3) 


where E' is the initial nucleon energy in the rest frame of the target. The differential 
spectrum is defined (nonrelativistically) over the energy interval 0 < E < E' . Note that 
the expression (3.5.1) reduces to the usual result for low-energy scattering: 


f(E,E’)*l/E' (3.5.4) 

The forward-to-backward scattering ratio is required for neutron scattering and is given by 
(ref. 46) 

F b (E') = 0.12 - 0.015£' + 0.4l/[l + e 4 ^'" 12 )] (3.5.5) 

where E in equation (3.5.5) is the laboratory energy (in GeV) before collision. 

The differential cross sections are normalized such that 


— =a(E')f(E,E') (3.5.6) 

where a{E') is the “appropriate” nucleon- nucleon total cross section. Obviously, we have 
neglected the inelastic processes that must yet be included so that a(E') in equation (3.5.6) 
is set equal to the total cross section to ensure particle conservation. The center of the mass 
angular distributions 9 cm is related to the energy change in the laboratory system by 

da E f da 

dn ~ AndE (3.5.7) 

(where Q denotes a solid angle) and is compared with the compilation of experimental data 
(ref. 49) in figures 32 and 33. These comparisons indicate that the present functions are 
reasonable. 


3.5.2. Nucleon- nucleus spectrum. The nucleon- nucleus differential cross section in 
Chew’s form of the impulse approximation (note that this is just the Born term of the optical 
model in ref. 50) is given by 


da 
dq 2 


- 2 bq *\F A {q 2 )\ 2 


« ce 


-2bq 2 e ~2a 2 q 2 /3 


(3.5.8) 


where b is the slope parameter of equation (3.5.2) averaged among nuclear constituents, q is 
the magnitude of momentum transfer, and a is the nuclear root- mean-square (rms) radius. The 
nuclear rms radius (ref. 50) in terms of the rms charge radius (in fermi) is given as 


a = 



(3.5.9) 
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where the rms charge radius (in fermi) is 


' 0.84 

(At = 1) 

2.17 

(At = 2) 

1.78 

(At = 3) 

1.63 

(At = 4) 

2.4 

(6 < A t < 

„ 0.82A < 1//3 + 0.58 

(A t > 16) 




> 


14) 


(3.5.10) 


The nuclear form factor is the Fourier transform of the nuclear matter distribution. Note that 
the above equation assumes that the nuclear-matter distribution is a Gaussian function. Such 
an approximation is reasonable for the light-mass nuclei but is less valid for At 'JS> 20. 

The energy transferred to the nucleus Et is restricted by kinematics to 


0 < E t < (1 -a)E' 


(3.5.11) 


where 

a = {At — l ) 2 / {At + l ) 2 

The energy-transfer spectrum is given as 


f(Et,E') 


4Atmc 2 (B + ^-) exp j^— 4Ttmc 2 (J3 + ^-)F'tj 
1 — exp — 4A<mc 2 (l — a)(B + ^p)!?^ 


Similarly, the nucleon energy after scattering E is restricted to 

aE’ <E< E' 


(3.5.12) 


(3.5.13) 


(3.5.14) 


The nucleon spectrum is given by 

AAtmc 2 {B 4- C-) exp [—4 Atmc 2 (B + ^p ){E' — E 1 )] 

f(E, E') j i 1 L (3-5.15) 

1 — exp — 4Ajmc 2 (l — a)(B + ^-)F' , J 

One should note that both equations (3.5.13) and (3.5.15) reduce to the usual isotropic 
scattering result at low incident energy. The differential spectrum is normalized as 


= a s (E’)f(E,E ') (3.5.16) 

where <r s {E') is the total scattering cross section obtained from equation (3.2.16). 

The results of equation (3.5.16) are compared with experiment (refs. 51 and 52) in figures 34 
to 37 (where <9 lab is the scattering angle in the laboratory). The comparison is rather good 
at the small angles when considering the simplicity of the present results. Also, shown in the 
figures are prior results from reference 2 showing considerable improvement in representing 
forward-scattered neutrons over the prior results. Much of the present discrepancy near 
forward scattering is due to errors in a s (E) to which the present spectra are normalized in 
equation (3.5.16). At broader angles, additional differences are due to the neglect of higher- 
order corrections to the impulse term. 

3.5.3. Nucleon nonelastic spectrum. The nonelastic differential cross sections (the 
inelastic process in which the nucleus is raised to an excited level that is ignored) use the results 
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of Bertini’s MECC-7 program (ref. 47). The nucleon multiplicities are given in tables 3.5.1 and 
3.5.2. We have required the multiplicities to be monotonic in energy, and thus the values m 
parentheses, which were obtained by scaling from lower and higher energies, are correct values 
and are used in the calculations. The results below 400 MeV were taken from Alsmiller et al 
(ref. 53), and the results for carbon, calcium, bromine, cesium, and holmium above 400 MeV 
are obtained by interpolation. The nonelastic spectra are represented as 




^N i _e^i Nq 

a< 1 - e~ E 'l a i E> [l + e -20(i-£/£')] 


(3.5.17) 


The first term of the summation represents the evaporation peak so that N\ (the number of 
evaporation nucleons) is taken from table 3.5.1 and the spectral parameter «i (in GeV) is taken 
from Ranft (ref. 54) as 


a lp 


■{ 


a ln 


(0.019 + 0.0017E')(1 - O.OOlAt) 
0.027(1 -0.00lA t ) 

(0.017 + 0.0017£')(1 - 0.00lA t ) 
0.023(1 - O.OOlAt) 


{E' <5GeV)l 
(E' >5 GeV) J 

(£'< 5GeVn 
(£' >5 GeV) J 


(3.5.18) 

(3.5.19) 


The second term is taken from Ranft (ref. 54) to represent the low-energy cascade particles as 


r 0.0035v^ (£' < 0- 1 GeV ) 1 

n 2p = | 0.007 v/Aj [0.5 + (1 + log 10 E') 2 } (0.1 < E' < 5 GeV) ► 

l 0.0245 s/A t {E' > 5 GeV) > 


n 2 n — 


/ 0.0042 s[A t 

< 0.007\/^[0.6 + 1.3(1 +log 10 £') 2 ] 

„ 0.032v^4f 


(E 1 <0.1 GeV) > 
(0.1 < E' < 5 GeV) > 
(£'> 5 GeV) . 


(3.5.20) 


(3.5.21) 


with the corresponding spectral parameters 


<*2p = 


«2n = 


( 0.11 + 0 . 01 E ')(1 - 0 . 001 A ,) 

0.16(1 - 0.001A t ) 

(O.l-t-O.01E')(l-0.OOlAt) 
0.15(1 - O.OOlAe) 


(£' < 5 GeV) | 
{E' > 5 GeV) J 

(£'< 5 GeV)) 
(£' > 5 GeV) J 


(3.5.22a) 


(3.5.22b) 


The third term in the summation is the balance of cascade particles after the inclusion of the 
quasi-elastic contribution. 

The quasi-elastic contribution is estimated by including the nuclear attenuation following 
the quasi-elastic event. The proton quasi-elastic cross section is 


a Q,pp = z t a PP + (-^t - Z*)' 

np 


a np | 


(3.5.23) 


OQpn = i A t ~ Z t)<7\ 
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and similarly for neutrons, 


^Q.nn — {At Z()a nn + Zt&np 

a Q,np ~ Zt(T n p 


The corresponding multiplicities are taken as 


(3.5.24) 


N Q,jk = e °'° 5 '^ J ‘ cr Q,jk / (3.5.25) 

’ £ 

where the exponential factor accounts for the attenuation of the quasi-elastic particles before 
they escape the nucleus. The balance of the cascade particles is contained in N% as 


N 3 = N c -N 2 - Nq (3.5.26) 

with an assumed spectral coefficient given by 

a 3 = a 2 / 0.7 (3.5.27) 

Results of the present formalism are shown in figures 38 to 51 in comparison with the 
calculations of Bertini et al. (ref. 46). Some further improvements in this parameterization 
need to be made. 


3.5.4- Light- fragment spectrum . The light-fragment yields per event are given in 
table 3.5.3 as obtained from Bertini’s MECC-7 calculations (ref. 47). These results are 
extrapolated and interpolated in energy and mass number. The corresponding mean energies 
are given in table 3.5.4. The mean energies are used in Ranft’s formula for nucleons and are 
similarly used for the light ions. 

3.5.5. Heavy- fragment spectrum. Following the direct ejection of nucleons in nuclear 
collision, the nucleus is left in a highly excited state that decays through particle emission. 
From a sudden approximation point of view, as proposed by Serber (ref. 41), the momentum 
distribution of the decay particles is governed by the fermi distribution prior to collision. The 
collective momentum of decay products and nuclear fragments is thus derived on the basis of 
combinatorial rules on the random ways in which a given fragment mass can be formed from 
the nucleon distributions prior to collision (refs. 55 and 56). The formulation of Goldhaber 
(ref. 56) is physically meaningful and simplistic. The momentum distribution is Gaussian in 
momentum space with a momentum width parameter given by 

a p = ct 0 [Af{A t - Af)/(A t - 1)] 1/2 

where ctq is the usual mean fermi momentum of the struck nucleus. However, the (Jq of nuclear 
fragmentation is found to be about 25 percent smaller than that observed in electron scattering 
experiments XreF. 39). The mean fermi momentum is a slowly varying function of nuclear mass. 

A slight modification of Goldhaber’s result is found to adequately represent the experimental 
results of Greiner et ah (ref. 39) given by 

CTp = 0.86[4<W20(.4f-l)] 1/2 
where the parameters b and 64 are given, respectively, by 

b = min (l!2/l^ 2 , 260) 


(3.5.29) 


(3.5.30) 


(3.5.28) 
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and 


^ = { A?- A, (Otherwise) } < 35M > 

A comparison of formulas (3.5.29) to (3.5.31) with experiments and the parameterization of 
Greiner et al. (ref. 39) is given in table 3.5.5. Clearly, the present formulas are quite accurate. 

The spectral distributions of the nuclear fragments in the rest frame of the struck nucleus 
prior to collision are given by 

^ = (dfyi7J £l/2 “'>(- £ / 2E ») < 35 - 32 > 

where cry is the fragmentation cross section and the energy parameter is 


Eq = 3<7p/2Af (3.5.33) 

The average energies E of various fragments obtained by equations (3.5.29) to (3.5.33) are 
compared with results of the Bertini model in table 3.5.6. Generally, the average energies 
predicted by the Bertini model are reasonably accurate, although some specific isotopes differ 
by a factor of 2 or more. 

3.5.6, Energy-transfer cross section. The energy-loss spectrum ipj(x : fi,£) of an ion 
fragment j (ref. 57) may be written as 


rpj{x,n,E) ss / 

J E 

where Aj is the fragment mass number, Cj( x ) is the fragment source, and E 1 is related to the 
distance to the boundary along the direction f 2 as given elsewhere (ref. 57). For distances far 
from the boundary, one may take Ey — oo. The cumulative energy-loss spectrum far from the 
boundary (Ey = oc) is 

Dj(x, E) — 47r f , ipj(x,rt y E r ) dE f (3.5.35) 

J E 

from which the distribution in linear energy transfer (LET) of energy deposit can be found. 
The total energy absorbed is given by 

D(x) = '£D j (x, 0) 

3 

EjOjp4> (3.5.36) 

3 

where Ej is the average energy of the fragment j, oj is the fragmentation cross section, 
p is the target density, and (j) is the effective nucleon flux initiating the fragmentation events. 
The energy- transfer cross section of the various fragment components is EjOj and is shown 
in table 3.5.7 for the Rudstam parameterization (present results), Bertini data (ref. 46), and 
experiments of the Heckman group (refs. 39 and 40) for comparison. Equations (3.5.34) to 
(3.5.36) also provide a basis for resolving the energy-transfer cross section into its various LET 
components. The LET components of equation (3.5.35) are shown in figure 52 for p = <f> = 1 
for all contributions with a fragment charge greater than 1. The two curves shown in the figure 
are for the Bertini data and the experiments of the Heckman group. Results obtained using 
our modified Rudstam formalism and the parameterized momentum distributions are virtually 
indistinguishable from the curve based on the LBL experiments. It is clear from the results 


/ \ 3 / 2 0 

f_ZLj 


(3.5.34) 
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shown in figure 52 that estimates of exposure from heavy ion recoil nuclei in tissue based on 
Bertini cross sections are generally low. 

4. Methods of Solution 

In an earlier paper (ref. 3), we proposed the use of a perturbation theory to develop 
a numerical method for solving the one-dimensional charged-particle transport problem. 
Although the resulting algorithm was greatly simplified compared with Monte Carlo algorithms, 
it still suffered from the need to manipulate large amounts of numerical data. In the present 
section, we show how the integral operators of that earlier work (ref. 3) may be numerically 
evaluated to eliminate the need to store and manipulate large amounts of numerical data, and 
at the same time we develop an algorithm that maintains a close relation to the physical field 
quantities. The resultant numerical solutions are compared with results obtained by analytical 
solutions for realistic interactions. 


4.1. Energy-Independent Proton Model 

The Boltzmann equation for proton transport in the straight ahead approximation is given 

fOO 

<t>(x,E) = f{E,E')<l>(x,E')dE' (4.1.1) 

Je 


l-W s{ - E) + ° 


where S(E) is the proton stopping power, a is the macroscopic interaction cross section which 
we presently take as energy independent, and f{E,E f ) is the production secondary-particle 
spectrum. Using the definitions 


r 

ip(x,r) 



S(E)<j>(x,E) 


(4.1.2) 

(4.1.3) 


and 

f(r,r') = S(E)f(E,E') (4.1.4) 


allows equation (4.1.1) to be rewritten as 

rx roc 

il’{x,r) = e~ ax 0(0, r + x) + / dze~ az / dr' f(r + z,r') ip(x — z,r') (4.1.5) 

J 0 Jr+z 


where the boundary condition is 


0(0 ,r) = S(E) 0(0, E) 


(4.1.6) 


A numerical algorithm for equation (4.1.5) is found by noting that 

fh roc 

ip(x + h,r) = e~ ah Tp(x,r + h) + J dze~ az J dr' f(r + z, r'z) ip(x + h - z, r' + z) (4.1.7) 
which can be simplified by using 

ip(x + h — z,r) ~ xp(x, r + h — z) + 0(h) (4.1.8) 

which yields 

ip(x + h, r) « e~ ah 0(x, r + h) + e~° h J dz J dr' f(r + z,r' + z) ip(x,r' + h) (4.1.9) 
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with the order of h 2 where h is the step size. Equation (4.1.9) is accurate for distances such 
that ah « 1 and may be used to relate the spectrum at some point x to the spectrum at 
x + h. Therefore, one may begin at the boundary (x = 0) and propagate the solution to any 
arbitrary interior point using equation (4.1.9). 

In the event that the boundary has a discrete spectrum such as 


then 


<t>{O i E) = 6{E-E 0 ) 


^(0 ,r) = S(r - ro) 


(4.1.10) 


(4.1.11) 


When discrete spectra are present at the boundary, the solution contains both singular and 
continuous components which we label ip s and ^ c , respectively. The corresponding singular 
term in the solution is then 

i>s(x , r) = e~ ax 6(r + x — r$) (4.1.12) 

whereas the continuous term satisfies 

tp c {x + h,r) = e~ ah xp c (x, r + h) 


rh roc __ 

4- J dze~ GZ J dr f f(r + z,r ; + z) [ips{% + h — z>r + z) 

+ V>c(^ + h — z,r f + z)] (4.1.13) 

The first term under the integral may be evaluated using equation (4,1.12) to obtain 

rh 

ip c (x + h,r) = e~ ah ip c (x, r + h) + e ~ a ( x+h ) / dz f(r + z, tq — x — h + z) 

Jo 

rh roc _ 

-f J dze~ GZ J dr f f(r + z,r f + 2 ) ^ c (x + h — z,r f 4 z) (4.1.14) 

The solution over small values of ( h — z) may be approximated as 

V> c (x + h — z, r) « e -<J (^“ 2r ) 7p c (x, r + h — z) - j- 0(h — z) (4.1.15) 

for which (see the appendix) 

ip c {x + h,r ) « e~ ah rp c (x, r + h) + e~ ax F^r, h,r 0 - x - ^)j 

+ e~ ah J dr'F(r,h,r')-ip c {x,r' + h) (4.1.16) 


where 


rh 

F(r, ft, r 7 ) = / dzf(r + z,r l ) 

Jo 

= F c \t{r + h),E , )-F c {E,E l ) 
where e(r) is the energy for range r and 

Fc(E,E')= f E f{E,E f )dE 
Jo 


(4.1.17) 


(4.1.18) 
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which is the cumulative secondary-particle spectrum. Note that equation (4.1.16) requires only 
one numerical integration per step in x. 

Discrete spectrum . Nucleon- nucleon scattering can be well- approximated (see 
section 3.5) by 

f{E, E') = ce- a{E '- E) (4.1.19) 

where c/a = a. This spectrum is related to the quasi-elastic spectrum of nucleon- nucleus 
reactions. Similar to this spectrum is 


f(r,r) = ce a ^ r (4.1.20) 

As a model problem, the spectrum of equation (4.1.20) is realistic and can be solved using 
perturbation theory. The first term is the uncollided beam term 

tpoix.r) = e~ ax b(r^ — r — x) (4.1.21) 

The first-generation term is 

^l(x, r) = xe~ ax ce~ a ( r ° ~ r ~ x ^ (4.1.22) 


and the higher-order terms are 

i>n{x,r) = ~x n e~ ax - - - -TT(r 0 - r - ar) n “ 1 e _a{r » _r_x) (4.1.23) 

n! (n — 1)! 

This problem is solved numerically using equation (4.1.16) and is compared with the analytic 
solution in table 4.1.1. The incident beam is for 500 MeV protons on a water shield with 
a = 0.01 cm 2 /g and a ^ 0.0123 cm 2 /g. 

As seen from the table, solutions with discrete spectra are limited in accuracy to 5 percent, 
independent of the depth of penetration. This error arises from the energy interpolation formula 
as the spectrum is highly discontinuous. Special interpolation methods could be developed to 
reduce this error greatly. 

4*1*2 . Continuous spectrum . For this test, a spectrum similar to a solar proton event 
is taken as 

^(0, r ) = e~^ T (4.1.24) 

The leading term in the perturbation theory is 

^o(x,r) = (4.1.25) 

with successive collision terms given by 

1 c 

ipn(x,r) = — x - — — —ip n ^i{x,r) (4.1.26) 

n (a + p) 

This problem is solved numerically and compared with the analytic result in table 4.1.2. It 
is seen from the table that the agreement for the two solutions in this case is generally within 
±1 percent. Clearly, high-quality numerical solutions are available for continuous spectra at 
the boundary. — 

The algorithm developed herein provides adequate solutions to proton beam problems 
with discrete spectral components and highly accurate solutions for typical space applications 
involving continuous spectra. The computation times for each of these test problems were less 


than several minutes on a CYBER mainframe, thus offering a very favorable comparison with 
Monte Carlo or previous methods based on the perturbation theory (ref. 3). 


4.2. Coupled Baryon Transport Methods 

The coupled baryon transport equations are of the form 

d 


d_ 

dx 


.,JLs(E) + <rjiE) } 4>j(x ,£) = £ [°° f jk (E,E')<t> k (x,E')dE' (4.2.1) 

Otl J T JO 


where vj is the range scaling parameter, S(E) is the proton stopping power, cr(E) is the total 
cross section, <f>j(x,E) is the differential flux spectrum of type j baryons, and fjk(E, E') is 
a differential energy cross section for redistribution of particle type and energy. Utilizing the 
definitions 


r = / dE r /S(E') 

J 0 

ipj{x,r) = S(E)<j>j{x,E) 


and 


f jk (r,r') = S(E)f jk (E,E') 

allows equation (4.2.1) to be rewritten as 


d d . 

S"iS + '' (r) 


=vf 


rpj(x, r) = 2J/ /jfc ( r - r ') 'Pkix, r') dr' 
k Jo 


which may be rewritten (refs. 3 and 4) as 

ipj(x, r ) = e~^> < r -*Vj(0, r + vjx) 


+ 


V / / e ^ ^fjkir + V]Z, t) ip k {x - t, r') dr'dz 

j? JO JO 

where the exponential is the integrating factor 

C jM) = J aj(r + Ujt')dt' 


(4.2.2) 

(4.2.3) 

(4.2.4) 

(4.2.5) 


(4.2.6) 


If the interactions are such that 


fjk( r . T ') = s jk 9( r ~ r ') (4.2.7) 

where g denotes the appropriate spectral function, then the solutions to equation (4.2.5) are of 
the form 

ipj(x,r) = x(x,r + Ujx) (4.2.8) 

To demonstrate how remarkable equation (4.2.8) is, we note that if x(arjr) is the solution to the 
neutron transport equation ( u n = 0), then xi x > r + I/ p x ) > s the solution to the proton transport 
problem independent of the functional form chosen for the stopping power. 

Rather simple numerical procedures follow from equation (4.2.6). Noting that the first-order 
nature of equation (4.2.1) allows ^(x, r) to be taken as a boundary condition for propagation 
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to larger x , one may approximate equation (4.2.6) as 


4’j{x + h } r) = e ^ r + Vjti) 

+ f f e-^{ T ' z }~fjk(r + VjZ,T f )4)(x + h — z,r , )dzdr r (4.2.9) 
T JO J 0 ' 


which may be used to develop a numerical stepping procedure. Equation (4.2.9) has provided 
the basis for a number of new transport codes for baryons of mass number greater than or 
equal to 1 (refs. 2, 4, and 5). These codes are now being extended to couple with the meson 
fields and to the negative baryon number fields. 

If h is sufficiently small such that 


<Tj(r) h << 1 (4.2.10) 

then, according to perturbation theory (ref. 3), 

4) k {x + h — z, r') « e“^ r ' + v k (h - z)] (4.2.11) 

which may be used to approximate the above integral of equation (4.2.9). 

For many cases of practical interest (e.g., accelerator studies), monoenergetic particle beams 
are used, and separation of the singular terms from the solution becomes convenient. The initial 
beam of type J particles of energy Eq (where tq = R(E$)) is taken as 


4’j(0,r) = 6jj6(r 0 - r) 

(4.2.12) 


and the solution is written as 



ipj{x, r ) = tpj 0 {x,r) + ipj{x,r) 

(4.2.13) 

| 

The corresponding singular terms are 



tf’jtote. r ) = e - *’* (r ’ x) <5(r 0 - r - v k x ) 6 kj 

(4.2.14) 


The regular terms of equation (4.2.9) for k — p may be written as 


K 

4’ p (x + h, r) = e ^ ( r ’^il>p(x, r + h) 


% 

+ f c/ze“^( r,z ) ^ [ f p j(r + z,r')[xl> j0 {x + h - z,r') 
JQ • Jr-^z 



+ xftj (x + h — z, /)] dr* 

(4.2.15) 

5 

and the regular terms for k = n are 



4> n (x + h, r) = e~ a " ip n {x, r) 


1 

+ Jo cZ2e_<7 " (r)2 E/ fnj(r,r’)[4>j 0 (x + h - z,r’) 


1 

+ ipj{x + h — z, r 7 )] dr 

(4.2.16) 



The singular contributions under the integrals of equations (4.2.15) and (4.2.16) can be 
evaluated with equation (4.2.14), and the approximations in equations (4.2.10) and (4.2.11) 
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can be applied to find 


i> p {x + h, r ) = exp[-cr p (r)/i] tp p (x, r + h ) 

4- exp | ~[(Jp{r) + <T p (r ' 0 )\ ^ | F pp (h , r, rg) <5 pj exp[-Cp(^o^)l 


+ 


+ 


+ 


and 


exp I — [<r p (r) + cr„(r 0 )]^ | F pn {h, r, r 0 ) 6 nj exp[-<T„(r 0 ) x ] 

J exp |— cr p (r) + a p (^r' +~^^F pp (h,r,r' + ^j ^ p (x,r'+h)dr' 

expj- <r p (r) + a n \r' + ^J^^F pn (h,r,r r ip n (x,r'+h)dr' 

(4.2.17) 


ip n {x + h,r) = exp[— <r n (r)/i] il> n (x,r) 

+ h f np {r, r' 0 ) exp j - [<x n (r) + a p (r ' 0 )] ^ | 6 pj exp[Cp(rg, x)] 

+ hf nn (r , r 0 ) exp |-[cr n (r) + <J n (r 0 )]^| 6 nj - exp[-<r„(r 0 )x] 

+ hj exp {~[°"( r ) + cr p( r, )]^| fnp(r,r')xp p dr> 

+ hj expj-[<r n (r) + <7„(r')]^| fnn{r,r') xp n (x,r') dr 1 (4.2.18) 
where t'q = rg — x — ^ and F is related to the cumulative spectrum F as given by 


Fij{h, r, 



fij{r + z,r')dz 


= Fij(r + h,r') - F^r,/) 


(4.2.19) 


with 

/•e(r) 

F ij( r i r 0 = y o fij{E, E f ) dE (4.2.20) 

e(r) is the energy associated with residual range r, and £' = e(r '). Equations (4.2.17) and 
(4.2.18) are evaluated by establishing an x-grid at which ipj(x m ,r) is evaluated where h is the 
distance between each successive evaluation. The integral over r f is accomplished by establishing 
an r-grid (and the corresponding E-grid) and using 


/ oo 00 f T t+l 

9 (r n ,r') tl>j(x m ,r')dr' ss Y^g n (r n ,f ( ) / ^(i m , r) dr' (4.2.21) 

n €=t» ,/r ' 

where 77 — (77 + r^ +1 )/2 and the series terminates at the highest E value in the r-grid. 
There is a spatially dependent discontinuity in the proton flux spectrum which requires right- 
and left-hand interpolation and integration. These discontinuities have been treated in the 
computational procedures. 
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4.3. Neutron Source 

The neutron transport equation in three dimensions is 


rOO 

[n • v + *„(£)] <t> n (x,n,E) = Y, J e\ n, n')^-(x, n', e') <m' dE' ( 4.3 


.3.1) 


Although the straight ahead approximation is adequate for most proton calculations, the 
neutron fields are more strongly affected by nonforward scattering components, particularly 
the low-energy neutrons. The reason is that the lower-energy neutrons have a greater range 
than the lower-energy protons because of the electric charge difference. Thus, a first-order 
correction to the straight ahead approximation may be applied by substituting the proton 
coupling in equation (4.3.1) by the straight ahead solution for 0 p (x,fi,JE7) so that 

roc 

Cn(x> n,E)= fn P (E , E\ n, n x ) 0„(x, E f ) dE 1 (4.3.2) 

J E 


The corresponding neutron transport equation is 


[n 


roc 

V+CTn(£)] <^n(x, n,E)= / fnn(E, E\ Q, O') O', E') dfi' dE 1 E) (4.3.3) 

JE 


The neutron source integral is treated in a fashion similar to that of equation (4.2.21). 


4.4. Target-Fragment Secondary Flux 

The target fragmentations produced in nuclear collision with the nucleon field must now be 
treated. The spectral parameters of the composite fragments are relatively independent of the 
projectile charge, energy, or direction. This leads to some simplifying assumptions so that 

% (x,n,E) = spy£,(x) J^fj(E')dE' ( 4 . 4 . 1 ) 

E-y = Rj l [Rj{E) + d(Q)} (4.4.2) 

where d(Q) is the distance from the boundary (ref. 57). The source of ions of type j is evaluated 
as ,00 

tow = 1 ^ n '- E "> dn> dE> ( 4 - 4 - 3) 

« Jo 

where aji(E') represents the fragmentation cross sections. The fj{E’) represents the spectral 
contributions averaged over all the target atomic constituents. In the present code, the distance 
to the boundary is assumed to be large. One could treat not only the boundary effects but the 
interface effects as well. 


5. Results 

As an initial validation of the present code, comparisons are made both with prior 
calculations using Monte Carlo methods and with experimental data. Fully three-dimensional 
Monte Carlo calculations have been made with the Bertini code as the nuclear cross section set 
augmented with low-energy neutron data. (See refs. 58 and 59 for a detailed discussion.) Energy 
absorption in a tissue slab for normally incident neutrons of energies 0.5, 2, and 10 MeV is shown 
in table 5.1. Also shown are the results of the present code. The results appear remarkably good 
when considering the crudeness of the straight ahead approximation for low-energy neutrons 
and the limitations on the present data base. Results for higher-energy neutrons are shown 
in figures 53 to 57. In each case, reasonable agreement with the results of Zerby and Kinney 
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(ref. 59) are obtained. Similar results are found for energetic protons as shown in figures 58 
to 61. The present data base changes the results in figures 53 to 61 by only a few percent 
when compared with the results in reference 2. In the present calculation, the first-generation 
proton spectrum is discontinuous for monoenergetic beams and is best handled by taking many 
energy points in the spectrum. However, the calculation time then becomes excessive. The 
present results were calculated using only 30 energy points. This process is adequate for space 
radiation, as shown in reference 60, but it is marginal for the present monoenergetic results. 
The use of numerical benchmark problems will allow us to understand the numerical procedures 
better. Such a benchmark has already provided some insight (refs. 60 to 62). 

The code has been used to calculate the dose behind various shields for typical space 
radiation. Three major solar-particle events of solar cycles 19 and 20 are represented in 
figure 62. The spectra as given in reference 63 have been used. The modification of the 
solar-event spectra at various depths in a lunar soil model is shown in figures 63 to 65. The 
importance of the buildup of secondary neutrons is clearly apparent in the February 1956 and 
November 1960 events and does not appear at all in the August 1972 event. The neutrons of the 
February 1956 event reach a stationary value between 25 and 100 g/cm 2 , as has been observed 
in our earlier calculations (ref. 37). The resulting dose within a 5-cm sphere of tissue- equivalent 
material is shown as a function of soil thickness in figure 66. The dose reduces only slowly for 
increasing the thickness beyond 20 cm. 

Galactic protons and their secondary neutron spectra behind varying thicknesses of alu- 
minum are shown in figures 67 and 68. The incident proton spectrum is that for a solar 
maximum according to the model of Adams et al. (ref. 64). It is clear that the neutron flux 
approaches a maximum near 50 g/cm 2 which is similar to the lunar soil results. Results for pen- 
etration of the martian atmosphere are indicated in figure 69. The potential use of polyethylene 
for controlling the neutron flux levels is indicated in figure 70. There are important geometric 
factors to be applied to all these results for which some detail is given elsewhere (ref. 65). 

6. Concluding Remarks 

The emphasis of the present code is on high-energy baryon transport, but such a code must 
adequately represent the low-energy neutrons in a reasonable way. It is seen from the present 
results that this representation has been accomplished in the present code. The calculation of 
100 to 400 MeV neutrons and protons on tissue is in reasonable agreement with a more complete 
Monte Carlo code. The primary advantage of the present code is computer efficiency while 
maintaining adequate accuracy. Future work will concentrate on improving the representation 
of the quasi-elastic peak, the low-energy neutron transport algorithm, and adding the effects 
of meson production to improve the comparisons further. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
December 21, 1988 
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7. Appendix 

Numerical Procedures 

In this appendix we consider the question of the appropriate numerical procedure for 
evaluating equation (4.1.5). The equation being solved is 


rx roc 

ip(x,r) = e~ <TX ^(0,r x) + / dze~ az I dr f(r + z, r f ) ift(x — z, r f ) 

Jo Jr+z 

which may be solved using perturbation theory for a monoenergetic beam as 

r ) = e ~ ax ${ r + x ~ r o) 

tpl(x,r) = e~ <7X [ dzf(r + z,rQ-x + z ) 

Jo 


(Al) 

(A2) 

(A3) 


with higher-order terms being obtained by repeated substitution into equation (Al). Note that 
if /(r, r f ) is a function only of (r — r 7 ), then 


^l(x, r) — xe GX f(r, r$ — x) 


(A4) 


Equation (A4) would hold for the quasi-elastic peak distribution for which 

/(r, r 7 ) « ce~ a ( r ~ r ) (A5) 

where a and c are constants. Equation (A4) does not follow for the quasi-elastic recoil particles 
for which 

/(r, r') « c'e~ ar (A6) 

It follows that equation (A3) can be written as 

ipl(x,r) « e~ ax [F(r + x } r 0 - x - Q) - F(r, r$ - x - Q )] (A7) 

where the choice for Q is not entirely clear but in some way represents the average z dependence 
of equation (A3) on the interval 0 to x. If c 7 and a are strictly constants, then equation (A7) 
is independent of the choice of Q. In general, the spectral function f(r,r r ) contains terms 
like those in equations (A5) and (A6) simultaneously so that whatever numerical solution is 
implemented, the character of both solutions (A4) and (A7) must be retained. 

The numerical solution to equation (Al) is rendered as 

ijj s (x -f ft,r) = e~ ah '4) s {x,r + ft) (A8) 


for the singular part and as 


ip c (x + ft, r) = e ah | r + h) + J dz F(r + z, ro - x - Q)e 

rh ro c 

+ J dz J F(r + z,r f + z)ip s (x>r / + h)dr f 


(A9) 


for the continuous spectral components. The first two terms of equation (A9) correspond to 
ij’(x,r) of the perturbation theory, and the requirement 


ipl(x + ft, r) = e 


-ah 


fa 


fh _ 

,(x,r + h)+ / dz F(r — z,tq — x — Q) i 
JO 


(A10) 
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i 


in 


must be met by any numerical procedure. One may show equation (A10) to be an identity for 
Q = h — z. We now inquire as to a suitable choice for Q. Clearly, computational accuracy of 
equation (A10) depends entirely on exactness to which the integral in (A10) is evaluated. 

We first note that the integral in equation (A 10) for the quasi-elastic peak has the value 

f h - 

/ F(r + z,tq - x — h — z) dz = hcexp[-a(rQ — x — h — r)] (All) 

J o 

for the exact value of Q. Assuming Q to be some fixed value results in 


th_ 

J F(r + z, ro - x — Q) dz = — exp[— a(ro — x - h — r)][e a< ^ - e a ^~ h ^] 


= hcexp[— a(ro — x — h — r)] 


1 + aQ — -ah + 0 (h 2 ) 
z 


(A12) 


Clearly, the error of equation (A12) is minimized by taking Q - h/2. Note that this value of Q 
would also be chosen on intuition, since it represents the values at the midpoint of the interval. 

The above question was investigated using the numerical calculations and an analytic 
solution for the quasi-elastic peak form of the secondary spectrum. The analytic solution 
is graphically presented in figure Al. The numerical solutions for Q = 0 and Q = h are shown, 
respectively, in figures A2 and A3. The corresponding errors are shown in figures A4 and A5. 
In each case, the errors mainly occur near the upper energy limit of the spectrum at each 
depth x. The Q - 0 solution is a slight overestimate of the flux at the highest energies and 
Q — h is a slight underestimate. In accordance with the result of equation (A12), we expect 
the errors of figures A4 and A5 to nearly cancel if the value Q = h/2 is used. The solution for 
Q = h/2 is shown in figure A6 with the corresponding errors in figures A7 and A8. Clearly, 
an adequate approximation is obtained using Q - h/2 , although it is clear from the present 
analysis that even greater improvements can be made. 
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Table 3.4.1. Present Correction Factors for Rudstam’s Formula 


AA 

Correction factor for — 

12 C 

16 0 

1 

1.3 

1.5 

2 

.5 

1.0 

3 

.3 

1.0 

4 

.1 

1.0 

5 

1.0 

1.5 

6 

.35 

.5 

7 


.5 

8 


.1 

9 


2.5 

10 


1.0 


Table 3.4.2. Comparison of Oxygen Fragmentation Cross Sections a of Reference 46 
With Experiments of Reference 40 and Parametric Results of References 42 and 43 


a f 

Fragmentation cross sections, cr, mb, from — 

Bertini 
(ref. 46) 

LBL 
(ref. 40) 

Rudstam 
(ref. 42) 

NRL 
(ref. 43) 

16 

7.0 

0.02 

8.7 


15 

85.1 

61.5 

61.0 

59.4 

14 

39.0 

35.4 

32.6 

32.2 

13 

13.9 

22.8 

29.7 

17.7 

12 

28.1 

34.1 

27.9 

36.0 

11 

5.0 

26.4 

31.4 

19.9 

10 

9.1 

12.7 

12.0 

11.0 

9 

1.0 

5.2 

7.1 

12.1 

8 

.2 

1.23 

2.1 

14.7 

7 

1.1 

22.2 

27.8 

19.4 

6 

3.8 

13.9 

18.0 

16.7 

Total 

193.3 

235.5 

258.3 

239.1 





Table 3.4.3. Comparison of Carbon Fragmentation Cross Sections a Measured in Experiments of 
Reference 40 With Two Simple Parameterizations 


(a) 12 C at 1000A MeV 


Fragmental 

tion cross sections, cr, mb, from — 

LBL 

Rudstam 

NRL 

(ref. 40) 

(ref. 42) 

(ref. 43) 

0.1 

6.7 

0 

55.3 

63.2 

69.0 

22.7 

28.0 

22.0 

5.8 

10.0 

15.2 

1.4 

4.8 

26.0 

18.9 

21.7 

20.7 

12.4 

14.7 

16.9 

116.6 

149.1 

169.8 


(b) 12 C at 2000A MeV 


Fragmentation cross sections, a, mb, from- 


L 

(ref. 40) 
0.09 
57.0 

22.7 
6.20 
1.6 

20.49 

14.8 


Rudstam 


(ref. 42) 

6.2 

60.4 

27.8 

10.4 

5.2 

24.4 

17.2 

151.6 


NRL 
(ref. 43) 
0 

58.5 

20.5 

14.2 
24.1 

19.9 
16.7 
































Table 3.5.3. Evaporated Ion Yields From Nucleon-Nucleus Collisions 
[Values in parentheses are for proton reactions] 



Ion yields at — 


500 MeV 

1000 MeV 

2000 MeV 

3000 MeV 

A t = 16: 





d 

0.111 (0.094) 

0.199 (0.237) 

0.257 (0.265) 

0.304 (0.311) 

t 

0.022 (0.029) 

0.024 (0.025) 

0.033 (0.025) 

0.029 (0.029) 

he 

0.018 (0.034) 

0.035 (0.043) 

0.037 (0.052) 

0.037 (0.048) 

a 

0.664 (0.400) 

0.720 (0.696) 

0.664 (0.624) 

0.640 (0.667) 

At — 27: 





d 

0.126 (0.130) 

0.245 (0.269) 

0.380 (0.396) 

0.442 (0.433) 

t 

0.028 (0.023) 

0.048 (0.052) 

0.063 (0.065) 

0.072 (0.069) 

he 

0.042 (0.035) 

0.067 (0.074) 

0.073 (0.091) 

0.083 (0.092) 

a 

0.370 (0.400) 

0.550 (0.566) 

0.597 (0.582) 

0.577 (0.577) 

A t = 65: 





d 

0.150 (0.171) 

0.379 (0.390) 

0.748 (0.766) 

0.935 (0.987) 

t 

0.031 (0.035) 

0.075 (0.068) 

0.145 (0.145) 

0.177 (0.191) 

he 

0.013 (0.014) 

0.039 (0.056) 

0.112 (0.124) 

0.166 (0.177) 

a 

0.124 (0.137) 

0.231 (0.231) 

0.373 (0.377) 

0.431 (0.441) 

A t = 100: 





d 

0.174 (0.183) 

0.456 (0.475) 

1.01 (1.02) 

1.44 (1.48) 

t 

0.028 (0.029) 

0.080 (0.081) 

0.207 (0.192) 

0.269 (0.273) 

he 

0.012 (0.017) 

0.055 (0.060) 

0.162 (0.185) 

0.249 (0.262) 

a 

0.158 (0.156) 

0.320 (0.339) 

0.490 (0.467) 

0.549 (0.540) 

A t = 207: 





d 

0.131 (0.152) 

0.536 (0.565) 

1.51 (1.57) 

2.54 (2.54) 

t 

0.038 (0.037) 

0.152 (0.163) 

0.415 (0.424) 

0.641 (0.644) 

he 

0.001 (0.002) 

0.017 (0.017) 

0.112 (0.106) 

0.211 (0.239) 

a . 

0.053 (0.063) 

0.195 (0.210) 

0.527 (0.514) 

0.751 (0.746) 
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Table 3.5.4. Mean Energies of Light Nuclear Fragments Produced in Nucleon-Nucleus Collisions 

[Values in parentheses are for proton reactions] 

Mean energies at— 



500 MeV 

1000 MeV 

2000 MeV 

3000 MeV 

5.55 

(6.19) 

7.91 

(7.89) 

9.55 

(9.81) 

11.1 

(9.80) 

6.10 

(6.40) 

8.33 

(8.69) 

9.71 (10.2) 

10.3 

(11.2) 

8.53 

(7.64) 

12.2 

(10.7) 

14.9 

(14.8) 

16.3 

(13.0) 

6.40 

(7.83) 

10.6 

(10.4) 

12.5 

(9.74) 

13.7 

(10.1) 

12.1 

(8.76) 

11.8 

(11.2) 

11.1 

(13.1) 

12.9 

(10.3) 

9.36 

(6.24) 

12.6 

(12.3) 

13.1 

(14.6) 

13.6 

(13.8) 

5.08 

(5.09) 

7.34 

(7.48) 

9.91 (10.5) 

11.6 

(12.0) 

6.87 

(6.90) 

8.61 

(8.92) 

11.1 

(11.9) 

13.5 

(13.7) 

9.57 

(9.42) 

10.8 

(11.2) 

14.3 

(14.8) 

17.2 

(17.4) 

9.16 

(9.54) 

10.8 

(11.1) 

13.0 

(13.9) 

16.6 

(13.7) 

10.5 

(10.8) 

12.5 

(12.8) 

13.4 

(14.1) 

14.4 

(14.5) 

12.7 

(13.4) 

13.2 

(13.6) 

13.8 

(13.8) 

14.5 

(14.6) 

4.24 

(4.32) 

5.67 

(5.70) 

7.92 

(7.91) 

9.67 

(9.58) 

8.25 

(8.30) 

9.66 

(9.76) 

12.1 

(12.3) 

14.4 

(14.2) 

9.88 (10.1) 

13.5 

(11.8) 

13.8 

(14.2) 

15.6 

(15.9) 

10.0 

(10.0) 

11.7 

(11.6) 

13.7 

(13.8) 

15.1 

(15.9) 

14.6 

(14.1) 

16.4 

(16.2) 

17.5 

(19.3) 

19.5 

(19.2) 

12.7 

(13.4) 

13.2 

(13.6) 

13.8 

(13.8) 

14.5 

(14.6) 

3.90 

(3.90) 

5.13 

(5.16) 

7.11 

(7.04) 

8.61 

(8.74) 

9.63 

(9.62) 

11.0 

(11.0) 

12.9 

(13.2) 

14.6 

(14.7) 

11.0 

(11.1) 

12.5 

(12.6) 

14.4 

(15.0) 

16.1 

(16.0) 

11.3 

(11.7) 

12.6 

(13.0) 

14.7 

(14.3) 

15.5 

(16.5) 

17.8 

(18.7) 

18.6 

(18.8) 

20.9 

(20.6) 

21.8 

(22.2) 

16.5 

(16.5) 

16.8 

(16.9) 

17.5 

(17.5) 

17.6 

(17.6) 

3.28 

(3.27) 

4.37 

(4.33) 

5.83 

(5.78) 

6.90 

(6.95) 

12.5 

(12.5) 

12.2 

(13.4) 

14.9 

(14.9) 

16.2 

(16.3) 

13.2 

(13.2) 

14.4 

(14.2) 

16.0 

(16.8) 

17.4 

(17.8) 

13.6 

(13.8) 

5.0 

(15.3) 

16.6 

(16.8) 

17.4 

(17.8) 

24.1 

(27.0) 

26.2 

(26.5) 

28.0 

(27.8) 

29.1 

(28.5) 

25.3 

(25.7) 

26.0 

(26.3) 

26.4 

(26.3) 

25.9 

(26.4) 



Table 3.5.5. a v for 16 0 Fragments Produced by 2.1 GeV Protons 


Value of (7 P , MeV/c, from 


Fragment 

Experiments 
(ref. 39) 

Present work 

Gremer 
(ref. 39) 

15 0 

94 ±3 

80.0 

83.8 

14 0 

99 ±6 

109.5 

113.1 

13 0 

143 ±14 

129.2 

133.5 

16 n 

54 ±11 

55.0 


15 n 

95 ±3 

80.0 

82.8 

14 N 

112 ±3 

109.5 

113.0 

13 n 

134 ±2 

129.2 

133.5 

12 n 

153 ±11 

143.4 

148.1 

15 C 

125 ± 19 

80.0 

82.8 

14 C 

125 ±3 

109.5 

113.10 

13 C 

130 ±3 

129.2 

133.5 

12 C 

120 ±4 

143.36 

148.09 

n c 

162 ±5 

153.45 

158.5 

io c 

190 ±9 

160.3 

165.6 

13 B 

166 ± 10 

129.2 

133.5 

12 B 

163 ± 8 

143.4 

148.1 

n B 

160 ±2 

153.5 

158.5 

io B 

175 ±7 

160.3 

165.6 

8 B 

175 ± 22 

165.5 

171.0 

n Be 

197 ±20 

153.5 

158.5 

10 Be 

159 ±6 

160.0 

165.0 

9 Be 

166 ±7 

164.24 

169.66 

7 Be 

166 ±2 

164.24 

169.66 

9 Li 

188 ±15 

164.24 

169.66 

8 Li 

170 ±13 

165.4 

171.0 

7 Li 

163 ±4 

164.24 

169.66 

6 Li 

141 ±7 

160.0 

169.66 

6 He 

167 ±20 

160.0 

165.0 






Table 3.5.7. Comparison of Fragment Energy- Transfer Cross Sections 
Ec 7 of Bertini With Experiments of Greiner/Lindstrom and Present Results 


A F 

Energy transfer cross sections, Ea , MeV-mb, from — 

Bertini 
(ref. 46) 

Greiner/Lindstrom 
(refs. 39 and 40) 

Present 

results 

16 

5.04 

0.0006 

0.26 

15 

60.6 

56.9 

56.4 

14 

48.8 

51.7 

47.6 

13 

37.6 

48.3 

62.9 

12 

85.8 

68.2 

55.8 

11 

37.9 

99.1 

117.9 

10 

52.8 

62.0 

58.6 

9 

6.5 

25.7 

35.1 

8 

2.5 

7.1 

12.1 

7 

6.11 

121.7 

152.4 

6 

31.4 

73.4 

95.1 

Total 

375.1 

614.1 

694.2 


Table 4.1.1. Ratio of Numerical Solution to Analytic Solution of Equation (4.1.23) 
for 500-MeV Protons on a Water Shield 


E , MeV 

Ratio for shield thickness, x, g/cm 2 , of — 

10 

20 

40 

60 

80 

0.1 

1.000 

1.000 

1.004 

0.998 

1.023 

19.5 

1.000 

.999 

1.004 

.997 

1.024 

120.9 

1.000 

.999 

1.002 

1.008 

1.046 

333.3 

1.003 

.994 

1.037 

.963 


454.1 

1.031 






Table 4.1.2. Ratio of Numerical Solution to Analytic Solution of Equation (4.1.26) 
for Continuous “Space” Proton Spectral Input on a Water Shield 



Ratio for shield thickness, x, g/cm 2 , of — 

E, MeV 

10 

20 

40 

60 

80 

0.1 

0.994 

0.995 

^ 0.998 

0.999 

0.999 

11.2 

.994 

.996 

.999 

1.001 

1.000 

36.3 

.997 

.998 

1.002 

1.003 

1.005 

118.1 

1.001 

1.003 

1.005 

1.006 

1.004 

383.9 

.997 

1.000 

.991 

1.000 

.996 





Table 5.1. Energy Deposition of 0.5-10 MeV Neutrons 


[Values in parentheses are from present calculations] 


Incident 
energy, MeV 

Depth, cm 

Energy deposition from Monte Carlo and present calculations, MeV 

Proton 

Heavy ion 

0.5 

0 1 

0.1107 


0.0073 



1-2 (1) 

.0986 

(0.0856) 

.0066 

(0.0076) 


4-5 

.0418 


.0028 



56 (5) 

.0331 

(.0135) 

.0018 

(.0026) 


9-10 

.0074 


.0003 



(10) 

10-11 

.0059 

(.0018) 

.0002 

(.0006) 


1415 

.0006 





(15) 

15-16 

.0007 

(.0013) 


(.0001) 


19-20 

.0002 





„„ (20) 


(.0002) 




20-21 

.0001 




2 

0 1 

0.2138 


0.0147 



1-2 (1) 

.1984 

(0.1887) 

.0133 

(0.0135) 


4-5 

.1539 


.0110 



(5) 

5-6 

.1349 

(.0818) 

.0105 

(.0092) 


9-10 

.0770 


.0054 



(10) 

10-11 

.0741 

(.0298) 

.0061 

(.0050) 


14-15 

.0301 


.0026 



(15) 

15-16 v 

.0240 

(.0226) 

.0019 

(.0024) 


19-20 

.0091 


.0008 



nn , ( 20 ) 


(.0097) 


(.0011) 


20 21 

.0114 


.0009 


10 

0 1 

0.3520 


0.0339 



1-2 

.3284 

(0.3377) 

.0342 

(0.0358) 


4 5 

.3220 


.0345 



5-6 (5) 

.2977 

(.2539) 

.0266 

(.0288) 


9 10 

.2674 


.0257 



(10) 

10-11 v 

.2661 

(.1759) 

.0225 

(.0214) 


14-15 

.2161 


.0203 



(15) 

15-16 v 1 

.2211 

(.1411) 

.0198 

(.0156) 


19-20 

.1635 


.0150 



(20) 

20-21 ' 

.1291 

(.1010) 

.0149 

(.0111) 
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BichseKref. 15) 

P a I m e r a n d Akhavan-Rezayat (ref. 20) 
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Northcliffe and Schilling(ref. 21) 



Figure 1. Calculated and experimental stopping powers in water for typical cosmic ray ions as a function of 
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Figure 2. The total proton-proton cross section of the present formalism compared with various experiments. 
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Fitrure 3. The total neutron-proton cross 


section of the present formalism compared with various experiments 






Lithium 


Figure 6. The total neutron-nucleus cross section according to the present formalism. 
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Figure 9. The total nucleon-copper cross section according to the present formalism, the theory of Townsend, 
Wilson, et al. (ref. 28), and various experiments. 
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Figure 10. The total nucleon-lead cross section according to the present formalism, the theory of Townsend 
Wilson, et al. (ref. 28), and various experiments. 






E, MeV 


Figure 12. The total neutron-nucleus cross sections of 8 0 and 13 A1 according to the prior data base (ref. 1) 
the present formalism, and an evaluated data base (ref. 27). 











E, GeV 


Figure 15. The neutron-carbon absorption cross section according to the present formalism compared with 
Letaw et al. (ref. 29) and various experiments. 
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Figure 16. The neutron-aluminum absorption cross section according to the present formalism compared with 
Letaw et al. (ref. 29) and various experiments. 




E, GeV 

Figure 17. The neutron-copper absorption cross section according to the present formalism compared with 
Letaw et al. (ref. 29) and various experiments. 
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Figure 18. The neutron-silver absorption cross section according to the present formalism compared with 
Letaw et al. (ref. 29) and various experiments. 



Figure 19. The neutron-lead absorption cross section according to the present formalism compared with 
Letaw et al. (ref. 29) and various experiments. 




Figure 20. The neutron-nucleus absorption cross sections at 1 MeV of the present formalism compared with 
various experiments. 



Figure 21. The neutron-nucleus absorption cross sections at 3 MeV of the present formalism compared with 
various experiments. 



Figure 22. The neutron-nucleus absorption cross sections at 5 MeV of the present formalism compared with 
various experiments. 



Figure 23. The neutron-nucleus absorption cross sections at 10 MeV of the present formalism compared with 
various experiments. 



Figure 24. The neutron-nucleus absorption cross sections at 14 NleV of the present formalism compared with 
various experiments. 



Figure 25. The neutron- nucleus absorption cross sections at 20 MeV of the present formalism compared with 
various experiments. 
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Figure 28. Oxygen fragmentation cross sections for baryon-14 isobars in proton collisions. 



Figure 29. Oxygen fragmentation cross sections for baryon-13 isobars in proton collisions. 





Figure 30. Oxygen fragmentation cross sections for baryon-9 isobars in proton collisions. 



Figure 31. Oxygen fragmentation cross sections for baryon-7 isobars in proton collisions. 






Figure 32. Neutron-proton differential elastic scattering cross section of present model compared with 
experiment. 



®cm 


Figure 33. Proton-proton elastic scattering differential cross sections of present model compared with 
experiment. 
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Figure 34. Neutron carbon scattering cross section according to present formalism, prior formalism 
and evaluated data (ref. 52). 
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Figure 35 Neutron scattering cross section for several elements according to present formalism, prior 
(ref. 1), and calculations of Fernbach (ref. 51) for 7 MeV neutrons. P 
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Figure 36. Neutron scattering cross section for several elements according to present formalism, prior formalism 
( ref . i), and calculations of Fernbach (ref. 51) for 14 MeV neutrons. 
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Figure 37. Neutron-uranium scattering cross 
and evaluated data (ref. 52). 
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Figure 38. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 100 MeV protons on oxygen. 



Figure 39. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 500 MeV protons on oxygen. 




Figure 40. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 1000 MeV protons on oxygen. 



Figure 41. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 3000 MeV protons on oxygen. 
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Figure 42. Nucleon cascade spectrum of Bertini compared with present work for neutrons produced by 
1000 MeV protons on oxygen. 
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Figure 43. Nucleon cascade spectrum of Bertini compared with present work for neutrons produced by 
3000 MeV protons on oxygen. 


Present 



Figure 44. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 1000 MeV protons on aluminum. 



Figure 45. Nucleon cascade spectrum of Bertini compared with cascade spectrum used to present work for 
protons produced by 3000 MeV protons on aluminum. 




Figure 46. Nucleon cascade spectrum of Bertini compared with present work for neutrons produced by 
1000 MeV protons on aluminum. 



Figure 47. Nucleon cascade spectrum of Bertini compared with present work for neutrons produced by 
3000 MeV protons on aluminum. 
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Present 



Figure 48. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 1000 MeV protons on lead. 



Figure 49. Nucleon cascade spectrum of Bertini compared with cascade spectrum used in present work for 
protons produced by 3000 MeV protons on lead. 



Figure 50. Nucleon cascade spectrum of Bertini compared with present work for neutrons produced by 
1000 MeV protons on lead. 



Figure 51. Nucleon cascade spectrum of Bertini compared with present work for neutrons produced by 
3000 MeV protons on lead. 
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Figure 52. Linear energy transfer (LET) cross section for fragmenting oxygen nucleus in water 



Figure 53. Dose in tissue to normal incident 60 MeV neutrons. 
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Figure 54. Dose in tissue to normal incident 100 MeV neutrons. 



Figure 55. Dose in tissue to normal incident 200 MeV neutrons 
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Figure 60. Dose in tissue to normal incident 300 MeV protons. 
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Figure 61. Dose in tissue to normal incident 400 MeV protons 






Figure 62. Time-integrated proton flux spectra for three anomalously large solar proton events. 
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Figure 65. Nucleon flux variation with energy for indicated /-values of regolith thickness for proton event of 
August 1972 (denoted by flare protons □ ). 
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Figure 66. Predicted 5-cm tissue dose equivalents for slab thicknesses between 0 and 100 cm in simulated lunar 




Figure 67. Galactic cosmic ray (GCR) proton flux variation with energy due to incident GCR protons for 
indicated t - values of depth in aluminum. 
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Figure 68. Neutron flux variation with energy due to incident GCR protons for incident t- values of depth in 
aluminum. 
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. Nucleon vertical flux at martian surface due to incident GCR protons. 



Figure 70. Attenuated surface neutron flux using a 2-g/cm 2 polyethylene shield. 
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